Appeared m Physica A, Vol. 233 (1996) pages 848-858 



Transport and Elastic Properties of Fractal Media 

Anthony P. Roberts 
Faculty of Environmental Sciences, 
Griffith University, Nathan, Qld. 4111, Australia 

Mark A. Knackstedt 
Department of Applied Mathematics, 
Research School of Physical Sciences and Engineering, 
Australian National University, Canberra ACT 0200, Australia 

October 1995 



Abstract 

We investigate the influence of fractal structure on material properties. We calculate the statistical 
correlation functions of fractal media defined by level-cut Gaussian random fields. This allows the 
modeling of both surface fractal and mass fractal materials. Variational bounds on the conductivity, 
diffusivity and elastic moduli of the materials are evaluated. We find that a fractally rough interface 
has a relatively strong influence on the properties of composites. In contrast a fractal volume (mass) 
has little effect on material properties. 

1 Introduction 

In the past decade a number of important natural and manufactured materials have been shown to exhibit 
microstructures with fractal characteristics. For example, it has been shown that sedimentary rocks have 
internal surfaces that exhibit fractal roughness and in some cases the fractal regime may include the 
entire pore volume [Q . This geometric complexity has a number of important implications for transport 
properties in porous rocks. In another example it was shown that fractured surfaces of metals are fractal 
in nature and proposed that the fractal dimensions of the surfaces could be correlated to the impact 
energy of the sample Q . Finally low- weight polymeric gel materials (aerogels) give a unique combination 
of properties; low-density and insulative, yet mechanically strong Aerogels are good examples of 
fractal materials, made up of clusters generated from the aggregation of primary particles [Q. 

While the concept of a fractal geometry provides a description of disordered material structure, it 
is not clear if the properties of these materials are affected by fractal characteristics nor how one may 
design and optimise fractal structure for a given application. Microstructural characteristics other than 
fractal dimension (e.g., connectivity, coordination) may be more important when deriving Structure- 
Property relationships. For example, in aerogels, the macroscopic connectivity at low density is the key 
microstructural characteristic associated with the mechanical strength of the material. 

To date there has been no explicit investigation of the influence of fractal structure on material 
properties. Theoretical studies |^, 0, H, || have relied on simplistic models of fractal microstructure. 
In general the relationship of such models to real materials is unclear, and the predictive ability of the 
models is qualitative in nature. In this paper we investigate the influence of fractal structure on material 
properties. We find that a fractally rough interface has a strong influence on material properties. A 
fractal volume, in contrast, has little effect on properties. 
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2 Evaluation of Bounds 



For general random composites there is no exact method of predicting properties. However a number of 
rigorous bounds, which depend on the microstructure of the composites, have been derived. These include 
bounds on the conductivity |l^ and the bulk pi and shear [|l2|, 13 moduli of composite materials 



(reviewed in Ref . [|14| ) . The bounds are expressed ||15| in terms of the volume fractions and properties of 
each of the phases and two microstructure parameters: 

^ r dr ds [\ ( P2ir)p2is) \ ... 

= 7^ — — duP2{u) [p3{r,s,t) (1) 

^PQ JoTjciSj-i \ P J 
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where ^3(7', s,t) is the three point correlation function [|l6j, = r"^ + — 2rsu and -P„(u) denotes the 
Legendre polynomial of order n. Until recently the use of the bounds has been restricted to 'particulate' 
media, such as uncorrelated overlapping spheres [l^ . Such models have limited utility in describing the 
microstructure of composite materials. We have shown that microstructure generated from level cuts of 
a random standing wave mimics the microstructure of a wide range of real composite materials including 
polymer blends [0, porous rocks |l^ and foamed solids ||l^. We have subsequently evaluated the 
bounds for a variety of level-cut microstructures and shown that it is possible to correlate quantitatively 



the effective physical properties of materials to their microstructure ||17|, . Berk 1 20 has shown how the 
level-cut scheme may be used to model fractal interfaces, and the extension to volume fractals is simple. 
Here we investigate the properties of these surface- and volume- fractal materials using rigorous bounds. 
This allows us to investigate the influence of fractal structure on material properties. 

A key quantity in the characterization of two phase materials is the two-point correlation function, 
P2if) ~ p(l — 1')7(^) + P^, where p is the volume fraction of the first phase and 7(r) is the normalized 
correlation function with 7(0) — 1, 7(00) = 0. Physically p2 gives the probability that two points a 
distance r apart will lie in the first phase. Two principle fractal types are evident in composites: surface 
and volume (or mass) fractals. Each can be characterized by the behaviour of the normalized correlation 
function "f{r). Surface fractal behaviour occurs when the interface between the phases is rough at all 
scales. In this case Bale and Schmidt have shown that 1 — j{r) ~ j-^-Ds ^ ^ q where 2 < D^. < 3 
is the surface-fractal dimension. Experimental techniques have determined that a variety of materials 
exhibit fractal surfaces with Dg > 2. These include crushed glass, zeolites 12^ , silica gels lignite 
coals 1^ , porous rocks Q , soils and composite steels |^ (see Korvin |25|| for a review) . On the other 
hand, volume (or mass) fractals occur when voids (or inclusions) are present at all scales. In this case 
7(r) ^ j.-i'i-D^) as r ^ cx) where < 3 is the volume-fractal dimension Fractal behavior of this 

type has been observed in silica gels [|[ ^ and soils . 



3 Model morphologies: The level cut Gaussian random field 

Random composites with a wide variety of morphological properties can be generated by taking level-cuts 
of a Gaussian random field (GRF). A simple definition of a GRF is 

yir)^J — '^cos{kik-r + (j>i) (3) 

where (j^i is a uniform deviate on [0, 27r) and ki is uniformly distributed on a unit sphere. The magnitude 
of the wave vectors ki are distributed on [0, 00) with a probability (spectral) density P{k) ( J P{k)dk = 1). 
A composite material can then be defined by taking the region in space where a < y{r) < (3 as the first 
phase, while the two regions contiguous to this y{r) < a; y{r) > (3 define a complementary second phase. 
In the case /3 = 00 a 1-cut material results; in the case (3 — —a a. 2-cut material. 
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The statistical correlation functions of these materials can be calculated [g9|, ^ 
fraction is given hy p = (27r)^i e~i* dt and the two point correlation function is 



The volume 
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Here g[r) = {y{0)y{r)) is the field-field correlation function and is related to the spectral density of the 
field 

' , sin /cr „ 

Pik)— dk. (5) 

kr 
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In terms of the normalized 2-point function j{r) — (p2('') ^ ) / {p ^ P^) it can be shown that ||2C 

1 — 7(7") ~ (1 — g{r))^ r — > & 7(r) ^ g{r) r —> 00. 



(6) 



Higher order correlation functions can also be evaluated for level-cut GRF media 1 31 , 32[ . The freedom 
in specifying a, (3 and P{k) (or g{r)) allows a wide variety of materials to be modelled. 



4 Surface fractals 

Recent microstructural studies of composite steels |^ have revealed a fractal boundary between a ferrite 
and pearlite phase. Using a scaling relationship between the area and perimeter of the pearlite grains it 
was found that df « 1.5. In three-dimensions this implies a surface fractal dimension ofDs — 1 + ~ 2.5. 
Small angle scattering studies have shown that lignite coal pi[ has Dg « 2.5 and sedimentary sandstones 
have 2.25 < Dg < 2.96 with many exhibiting Dg = 2.5 ± 0.1 [Q. Materials with such a fractal surface can 
be generated in the GRF scheme by using, for example, g{r) — e^*"/' where / is a pore/grain length scale 
which we normalize to unity. In this case 1 — 7(r) ^ r 2 as r — > so that Dg = 2.5. The corresponding 
spectral density is 

To account for the fact that physical surfaces only exhibit fractal scaling down to some finite length 
scale (e.g. the lattice constant) a cut-off length can be introduced in the model by setting P{k) = for 
k > K (and scaling so / P{k)dk = 1). This leaves the large scale microstructure unchanged and filters 
out the high frequency (rougher) sinusoids in Eqn. (||) above a length scale 2t:/K. The magnitude of 
K is estimated from experimental studies. In ferritic-pearlite steels [|| the minimum scale reported is 
.lfj,m and the grain size is O(50/xm) so K ^ 27r50/.l = 3.1 x 10^. In sandstones [Q the minimum scale 
measured is .5 x 10~^iim and the pore-size is around 10/im so K ^ 27rl0/.0005 = 1.3 x 10^. We assume 
K = 00 provides a reasonable model of fractal steel and sandstone interfaces. To clearly demonstrate 
the effect of a fractal surface on material properties we also consider materials with K — 8 (providing a 
smooth interface with Dg — 2). 

Cross sections of the surface fractal material are shown in Figs. |^ (a) & (b). The morphology of the 
2-cut field is very similar to that observed in steel composites The roughness of the interface for the 
1-cut case is clearly visible in Fig. |^. The spectral density and associated 2-point correlation function of 
each of the materials are shown in Fig. Ten two-dimensional realizations of each model are calculated 
from Eqn. (^) for the case a = 0& /3 = oo {p — The average value of 7(r) is represented by symbols 
in this figure: the fractal scaling of 7(r) is evident. 

To determine the effect of the fractal microstructure on the conductivity and elastic moduli of each 
material we evaluate the microstructure parameters for both models (K — 00,8). The results are given 
in Tables |l]& || and bounds on the shear moduli are shown in Fig. ^ for the case where the the shear and 
bulk moduli are both equal to 1 in phase 1 and in phase 2. Note that the lower bound vanishes at this 
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(c) 1-cut £>„ = 1.5 (d) 2-cut L>„ = 1.5 



Figure 1: Cross-sections of the fractal composites generated using Eqn. The black region corresponds 
to Phase 1 at a volume fraction p = 0.2. In the 1-cut case 0.84 < y{r) < oo and in the two cut case 
-0.25 < y{r) < 0.25. 



contrast. Curves for the bulk moduli and conductivity are qualitatively similar. The shear moduli of the 
fractal and smooth 1-cut models are quite similar. In contrast the fractal 2-cut model has a significantly 
reduced resistance to shear when compared to the smooth 2-cut model: a rough interface decreases the 
elastic strength of a material if the inclusion phase is stronger than the bulk phase. This result is similar 
to a trend observed in composite steels: an increase in Dg leads to a decrease in impact toughness. 
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Figure 2: The level cut media, (a) A 1-cut surface-fractal interface {Dg — 2.5) at p = 0.22 (0.75 < y{r) < 
oo); (b) A 2-cut volume-fractal {Dy = 1.5) at p = 0.2 (-0.25 < y{r) < 0.25). 



0.40 




Figure 3: The spectral density and normalized 2-point correlation function of the 1-cut surface fractal 
media (p = ^). The symbols represent computational measurements of j{r) averaged over 10 cross- 
sections. 
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5 Volume fractals 

Aerogels provide good examples of volume fractal materials with w 2 0] . The formation process has 



been modelled by the cluster-cluster aggregation model [g3|, 28| which gives Dy « 1.7 — 1.8. To see 
how the level-cut scheme can model volume-fractal materials consider j{r) w C{r/l)~^ for r — > oo. The 
volume of solid within a radius i? is Q 

V{R) = I dvp2{r)/p= / drAnr'^ {{1 - p)j{r) + p) 
Jo Jo 

Now the the former term dominates for R < Rs — ^ ^^(3!^^^ ^ I ^ p^^/^l where Rg is a saturation 

scale and V{R) ~ R^^^ = R^" . Thus the volume fractal dimension is _D„ = 3 ~ A. A simple field-field 
correlation function which gives rise to this behaviour in the level-cut GRF is 

+ . (9) 

This field-field correlation function actually leads to a spectral density 

„/, s 1 1 , X sin fcr , , , 

^(fc) - ^ AiTr^g{r)-^dr (10) 

which is not strictly positive. Therefore we define a new P{k) which is set to zero for all k beyond the 
first point K at which P{k) < (and re-scale to ensure / P{k)dk = 1). This modification does not 
change the behaviour of the model at large r. In aerogels / (which we normalize to unity) is related to 
the length scale of the monomers. In experiments the saturation scale is O(lOnm) Q and in a recent 
model of a silica gel based on the cluster-cluster aggregation algorithm the saturation scale is around 
six times the particle diameter at concentrations of c = 0.05 py]. This is consistent with Eqn. (N): 
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Table 1: The microstructure parameters the fractal and regular composites for 1-level cut GRF's ((3 — oo). 





Model 




Ds = 


= 2.5 


Ds = 


= 2.0 


= 


= 1.5 




= 3.0 


p 


Ci 


m 


Ci 


Vi 


Ci 


Vi 


Ci 


m 


0.05 


0.294 


0.247 


0.215 


0.165 


0.220 


0.172 


0.197 


0.148 


0.10 


0.319 


0.276 


0.258 


0.212 


0.263 


0.219 


0.243 


0.197 


0.15 


0.342 


0.305 


0.294 


0.253 


0.299 


0.260 


0.281 


0.239 


0.20 


0.366 


0.332 


0.327 


0.291 


0.330 


0.296 


0.316 


0.279 


0.25 


0.388 


0.360 


0.358 


0.327 


0.360 


0.331 


0.348 


0.317 


0.30 


0.411 


0.390 


0.387 


0.363 


0.388 


0.365 


0.380 


0.354 


0.35 


0.433 


0.415 


0.416 


0.396 


0.417 


0.398 


0.411 


0.391 


0.40 


0.456 


0.442 


0.444 


0.431 


0.444 


0.431 


0.441 


0.427 


0.50 


0.500 


0.500 


0.500 


0.500 


0.498 


0.496 


0.500 


0.500 



Table 2: The microstructure parameters the fractal and regular composites for 2-level cut GRF's {(3 — 
-a). 





Model 
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= 2.5 
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= 2.0 


D, = 


= 1.5 


D,. -- 


= 3.0 


p 


Ci 


m 


Ci 


'71 


Ci 


Vi 


Ci 


m 


0.05 


0.402 


0.351 


0.786 


0.613 


0.768 


0.589 


0.802 


0.627 


0.10 


0.415 


0.369 


0.706 


0.516 


0.689 


0.500 


0.727 


0.535 


0.15 


0.431 


0.387 


0.656 


0.480 


0.636 


0.450 


0.684 


0.501 


0.20 


0.449 


0.402 


0.628 


0.473 


0.615 


0.470 


0.657 


0.491 


0.25 


0.463 


0.430 


0.612 


0.478 


0.602 


0.473 


0.642 


0.495 


0.30 


0.481 


0.451 


0.605 


0.491 


0.596 


0.488 


0.636 


0.508 


0.35 


0.497 


0.474 


0.603 


0.510 


0.596 


0.506 


0.634 


0.526 


0.40 


0.515 


0.495 


0.607 


0.533 


0.601 


0.529 


0.636 


0.546 



Rs p~^l-^ — 0.05^^/^'^ — 5.8. Although the level-cut materials generated using Eqn. (H) has fractal 
scaling over the correct range of r, it does not necessarily follow that this is an appropriate model for 
aerogels. It does, however, allow the dependence of properties on to be gauged within the level-cut 
scheme. 

We choose a fractal A — 1.5 (Dy = 1.5) and non-fractal case yl = 3 (so the second term in Eqn. 
dominates and Dy — 3) to determine the effect of volume fractal behaviour on material properties. The 
functions P{k) are shown in Fig. ^. We have generated realizations of the materials using Eqn. (^). 
Cross-sections of the model for the case Dy = 1.5 are shown for the 1-cut and 2-cut cases in Figs. |l| (c) & 
(d). A three dimensional representation of the 2-cut interface is shown in Fig. ^ The 2-point function of 
each model has been measured for the case P ~ ^ and averaged over 50 cross-section realizations of each 
model. The results, plotted in Fig. ^, show that the level-cut GRF model has volume-fractal scaling. 

The microstructure parameters for both Dy — 1.5 and Dy = 3.0 are given in Tables |l| & || and bounds 
on the conductivity are shown in Fig. ^ for the case f7i^2 = 1,0 (the lower bounds vanish at this contrast). 
As we would expect the sheet-like nature of the structures in the 2-cut media are much better conductors 
than the node/bond-like structures present in the 1-cut media. We find virtually no difference in the 
fractal and non-fractal materials in both the 1-cut and 2-cut cases. This surprising result suggests that 
the properties of volume-fractal composites (such as aerogels) are not explicitly dependent on the fractal 
dimension. 
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1 2 3 4 5 0.1 1.0 10.0 

k r 

Figure 5: The spectral density and normalized 2-point correlation function of the 1-cut volume fractal 
media [p — i). The symbols represent computational measurements of 7(r) averaged over 50 cross- 
sections. 
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6 Conclusion 

We investigate the influence of fractal structure on material properties. We have calculated rigorous 
bounds on the conductivity and elasticity of fractal media generated using the level-cut random field 
model. The behaviour of the bounds indicates that a fractal interface plays a minimal role in the properties 
of 1-cut media. For the two-cut model, which mimics the microstructure of both foams |]l9| and porous 
rocks a much stronger influence is observed. In contrast, varying the volume-fractal dimension of 
both 1-cut and 2-cut media has little effect on the property bounds. The latter result indicates that 
the remarkable properties of aerogels are influenced more by the fact that they contain very well 
connected structures at high porosities, rather than their fractal characteristics |§, In the future 
we shall utilise a range of microstructural studies |^ to develop a more appropriate model of aerogel 
structure. This will allow a more rigorous comparison between model and experimental properties. 
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